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ABSTRACT 

Blazars are characterized by large amplitude and fast variability, indicating that the 
electron distribution is rapidly changing, often on time scales shorter than the light 
crossing time. The emitting region is sufficiently compact to let radiative losses dom- 
inate the cooling of high energy electrons. We study the time dependent behaviour 
of the electron distribution after episodic electron injection phases, and calculate the 
observed synchrotron and self Compton radiation spectra. Since photons produced 
in different part of the source have different travel times, the observed spectrum is 
produced by the electron distribution at different stages of evolution. Even a homo- 
geneous source then resembles an inhomogeneous one. Time delays between the light 
curves of fluxes at different frequencies are possible, as illustrated for the specific case 
of the BL Lac object Mkn 421. 
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1 INTRODUCTION 

Variability is one of the defining properties of blazars, char- 
acterized by variations of their flux even of two orders of 
magnitude in time scales of years, and smaller changes, but 
■ still up to factor 2, in hours/days time scales. Observa- 
tional efforts to characterize the variability behaviour, which 

i 

started soon after the discovery of blazars, have been inten- 
sified in recent years, to study very fast fluctuations and 
possible correlation of fluxes at different frequencies (for a 
review see Ulrich, Urry & Maraschi, 1997). This will hope- 
fully shed light on the location of the emitting regions and on 
the nature of the acceleration mechanisms. Theoretical ef- 
forts, however, lag behind observations: the classical paper 
on time dependent synchrotron and Compton flux is still 
Kardashev (1962), which solves the continuity equations for 
the particle distribution in cases where the cooling time t coo i 
can be considered longer than the light crossing time R/c 
for any energy. Recently, Mastichiadis & Kirk (1997) and 
Atoyan & Aharonian (1997), apply the results of solving the 
continuity equation for the electron distribution to the vari- 
ations seen in the BL Lac object Mkn 421 and the galactic 
superluminal source GRS 1915+105. In both cases, only the 
variability on time scales longer than R/c were studied. On 
the other hand, significant variability often happens on time 
scales extremely short in highly luminous (and compact) ob- 
jects. This implies that the cooling time for the highest en- 



ergy electron may well be shorter than R/c, even once the 
effect of Doppler boosting and blueshift is accounted for. 
This is also indicated by the behaviour of the light curves 
as seen, e.g., in the X-rays and in the optical, showing a 
quasi-symmetric behaviour, with rise and decay time scales 
approximately equal (see i.e. Urry et al. 1997, Ghisellini et 
al. 1997, Massaro et al. 1996, Giommi et al. 1998), indicating 
that both times are connected to the light travel time across 
the source R/c, and therefore suggesting that the cooling 
times of the emitting electrons are shorter. This in turn im- 
plies that the electron distribution, at least at these energies, 
is significantly changing on time scales shorter than R/c. 

In these cases the knowledge of the time evolution of the 
particle distribution of the emitting electrons is not enough 
to construct theoretical light curves for the observed fluxes, 
since light travel time effects play a crucial role. The observer 
in fact sees, at any given time, photons produced in different 
part of the source, characterized by a particle distribution 
of a different age. The observed flux is then the sum of the 
emission produced by different particle distributions. 

In §2 of this paper we present our method of solving 
the continuity equation for the emitting particle distribu- 
tion, together with the assumptions made. In §3 we discuss 
illustrative examples of the time behaviour of the particle 
and photons distributions, having care to evidentiate the ef- 
fects introduced by the different photon travel times. Besides 
considering simple illustrative cases where the injection of 
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'fresh' relativistic particles is assumed to occur simultane- 
ously throughout the source, we also study more realistic 
cases which simulate the injection occuring when a shock 
front travels down a region of a jet. In §4 we apply our 
model to Mkn 421, to show how it is possible, even in our 
simplified model, to explain the time lags observed in the 
hard and soft X-ray light curves. Finally, in §5 we discuss 
our findings. 



modified for our purposes, since we are dealing with a conti- 
nuity equation with injection and escape terms and no heat- 
ing. The Chang & Cooper (1970) scheme allows to find more 
stable, non-negative and particle number conserving solu- 
tions. Moreover that scheme significantly reduces the num- 
ber of meshpoints required to obtain accurate solutions, and 
it is also appropriate for including heating terms (Ghisellini, 
Guilbert & Svensson, 1988). 

For all our runs we use an energy grid with equal loga- 
ritmic resolution: the energy meshpoints are defined as 



2 THE MODEL 
2.1 Assumptions 

We assume that the emission is produced by a distribution of 
relativistic electrons injected in a region of typical dimension 
R embedded in a tangled magnetic field B, at a rate Q(y) 
[cm -3 s -1 ] (7 is the Lorentz factor). Electrons lose energy 
by emitting synchrotron and synchrotron self-Compton ra- 
diation (SSC); they can also escape from the emitting region 
in a time scale t esc , assumed to be independent of energy. 

Our main purpose is to apply our calculations to the 
study of the blazars short time scales variability, espe- 
cially in frequency bands around and above the peaks of 
blazars spectral energy distribution. In these spectral re- 
gions synchrotron self-absorption is negligible, and conse- 
quently we neglect this heating effect in the kinetic equa- 
tion (e.g. Ghisellini, Guilbert & Svensson, 1988). Photon- 
photon collisions, producing electron-positron pairs, are also 
neglected, since we will deal with sources of small compact- 
ness £ = La T /(Rmc 3 ) (Maraschi, Ghisellini & Celotti, 1992) 

The continuity equation governing the temporal evolu- 
tion of the electrons distribution N(j, t) [cm -3 ] is 



^ = » W7,*W7,*)] + 0(7,*)-^ 
at o-y t esc 

where 7 = 7s + 7c is the total cooling rate, given by 

4 arc 



7=o 



3 m e c 2 



[u B + u Tad [rt,t)]r 



(i) 



(2) 



where ot is the Thomson cross section, Ub is the magnetic 
field energy density and U rac t{'y,t) is the energy density of 
the radiation field. 

The cooling time scale is assumed to correspond to syn- 
chrotron and self-Compton radiative losses: 

_ 3m e c 2 

tc °° l - AtTTCTfpB + U rad ] (6 ' 

This relation can be rewritten by assuming that all radia- 
tion energy density is available for scattering, and that it 
corresponds to the compactness £ = 4tt RarUrad/ (m e c 2 ): 

tcool _ 37T . . 

Wc = 741 + U B /Urad] 1 ' 

Note that for compactness values in the range 10 _4 -10 _1 , 
typical for the IR to 7-ray emission of blazars, the high 
energy electrons cool faster than R/c. 



2.2 Numerical method 

We numerically solve equation ([!]) adopting the fully implicit 
difference scheme proposed by Chang & Cooper (1970), 



7j — 7 m * n 



(5) 



where j m ax is the meshpoints number, and the energy inter- 
vals are A"fj = 7j+i/2— 7^-1/2- Quantities with the subscript 
j ±1/2 are calculated at half grid points. In the performed 
simulations a grid of 200 points has been used both for par- 
ticle energy and photon frequency. 

In order to discretize the continuity equation we define 



nTi+l _ -i AT i + 1 
*j±l/2 - 7j±l/2-< V 3 -±l/2 

and equation (jl|) can be written as 

-1 Aji C^+l zri + l 
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In this specific case we have Nj + \/2 = Nj+i and Nj -1/2 = 
Nj, according to the prescriptions of Chang & Cooper 
(1970). We can now rewrite the continuity equation as 



V3jN}Xl + V2 J iV; +1 + Vl 3 N]t\ = S) 
where the V coefficients are 



Vlj = 



V2 3 - = 1 + — 1 + 



At At7,_j. 



j-1/2 



V3, = 



and 



tesc ^7j 

At7j+i/2 

Mi 



(9) 



(10) 



S) = Nj + Q) At . 



(11) 



The system of equations (|9|) forms a tridiagonal matrix, and 
it is solved numerically (e.g. Press et al. 1989). We tested 
our method with the analitic solutions given by Kardashev 
(1962) in the case of synchrotron radiation only and injection 
of a constant power-law and monoenergetic distributions. 

We then calculate the synchrotron emissivity e s (u, t) 
[erg s _1 cm -3 sterad -1 Hz -1 ] of each distribution N^^ti) 
with 



= ^ / dntN(n,t)P.{y,i) 



(12) 



In the above formula P s (v,^) [erg s _1 Hz -1 sterad -1 ] is 
the single particle synchrotron emissivity averaged over an 
isotropic distribution of pitch angles 



Ps (f,7) 



3\/3 otcUb 2 

7T Vb 



(13) 



{K 4/3 (t) K l/3 (t) ~ |t [Kl /3 (t) - K 2 1/3 (t)]} 
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(Crusius & Schlickeiser 1986; Ghisellini, Guilbert & Svens- 
son 1988), where t = u/(3j 2 ub), vb = eB/(2Tvm e c), K a {t) 
is the modified Bessel function of order a. We then calcu- 
late the synchrotron radiation field I s (u,t) [erg s _1 cm -2 
sterad -1 Hz -1 ] using the transfer equation 



k(u,t) 



-k(v,t)Rl 



(14) 



here k(u,t) [cm 1 ] is the absorption coefficient (e.g. Ghis- 
ellini & Svensson, 1991) 



k(u,t) 



8ivm e u 2 



«- JV( 7 ,t) d 
7P d'y 



[ypP(-y,u) 



(15) 



where p — (j 2 — l) 1 ' 2 is the particle momentum in units 
of m e c. In order to calculate the inverse Compton radiation 
field, we make the assumption that the synchrotron radia- 
tion insta ntan eously fills the whole of the emitting region 
(see also §2.3) and we take into account the Klein- Nishina 
decline using the following approximation (Zdziarsky, 1986) 



ut for 7a; < 3/4 
for jx > 3/4 



(16) 



where x = hu/m e c 2 . We approximate the energy density of 
the synchrotron radiation U ra d, sy n{y,t) by setting 



4-7T 

Urad,syn (7; ^) — 

C 



47) 



dv Is (u, i) 



(17) 



and the integration limits are defined as u m i n = u Sl min 
(minumum synchrotron emitted frequency) and Umaxi^y) = 
min [vs.max, 3m e c 2 /(4/i7)] . 

The inverse Compton emissivity can be calculated, with the 
above assumptions, following Rybicki & Lightman (1979) 



e c (ui) 



(JT 

T 



duo 

UQ 



" y iV( 7 )/(^i)-I s M(18) 



u 



where uo is the frequency of the incident photons, u\ is the 
frequency of the scattered photons, [3 — v/c, and the inte- 
gration limits are 



71 = max 



72 



f— V 



1 ^min 



^max 1 



3 m e c 



4 hvo 



(19) 



(20) 



and where u™ ax and Uq 11 " axe the extreme frequencies of the 
synchrotron spectrum. The function f(uo,u\) is the spec- 
trum produced by the single electron, scattering monochro- 
matic photons of frequency uo (see e.g. Rybicki & Lightman, 
1979) 

(1+0)% -(1-/3); ^f<^<l 

f( Va , VX )=\ (l+0)-%(l-P); 1<^<1±| (21) 







otherwise 



The inverse Compton radiation field is simply obtained 
as Ia{ui) = e c (n) R- 



2.3 Slices 

The knowledge of AT (7, t) is not sufficient to properly re- 
produce the variability behaviour of blazars: observations in 
the optical band and at higher frequencies strongly suggest 
that the injection and/or the cooling times can be shorter 
than the light crossing time R/c. In this case the particle 
distribution evolves more rapidly than R/c and the observer 
will see a convolution of different spectra, each produced in 
a different part of the source. Initially the observer only sees 
the emission coming from fresh electrons located in the re- 
gion ('slice') closest to her/him; then the innner parts of the 
source become visible, also showing 'young' spectra, while 
electrons in the front slices are evolving. After a time R/c 
all the emitting region will be visible: the back of it with 
fresh electrons and the front of it with older electrons. 

In order to reproduce this effect we can divide the source 
of size R into n slices of equal thickness R„i and then sum 
the contribution of each slice to obtain the correct observed 
spectrum. Note that the number of slices has to be chosen in 
order to have R„i /c < tmin, where t m i n is the shorter among 
the relevant time scales (cooling, injection and escape time): 
in this way each slice can be considered an homogeneous re- 
gion, while the different slices are observed in a different 
evolution time, resembling an inhomogeneous emitting re- 
gion. 

We consider the case in which all slices have equal vol- 
umes, and assume a 'cubic' geometry, with the line of sight 
placed at 90° with respect to one face of the cube. This 
angle, besides corresponding to the simplest case, is trans- 
formed, in the lab frame, into a viewing angle of ~ 1/r, 
if the source moves with a Lorentz factor F, which is ap- 
propriate for blazars. Extension to different geometries (i.e. 
cylinder, sphere, etc) is trivial, by properly weighting each 
slice volume. We will see that even in the simplest, 'cubic', 
case, time-lags among light curves at different frequencies 
are clearly observable. 

We note that the light crossing time effects should be 
taken into account also for the correct calculation of the 
inverse Compton radiation: synchrotron photons in fact re- 
quire a time R/c to fill the whole of the emitting region. 
This effect, relatively important for transient high energy 
phenomena in Compton dominated sources, is however neg- 
ligible for magnetic dominated objects, and we neglect it in 
the present paper, due to the computing time it needs to be 
properly calculated. We plan to improve our numerical code 
in the future. We expect that this effect will cause a time 
delay between synchrotron and self-Compton fluxes. 



3 SIMULATIONS 

3.1 Gaussian injection with ti„j <C R/c 

We first examine the case of a narrow Gaussian electron dis- 
tribution continuously injected for a time ti n j C -fi/c in a 
region of size R — 10 16 cm embedded in a tangled mag- 
netic field B = 1 Gauss. The injected distribution peaks at 
7 = 10 5 . We measure the total injected power with the corre- 



sponding compactness, defined as 



in this case we use 



t 



10" 



l in j = Li„ja T /(Rm e c ): 
The injection stops at 



i?/(10c), and we set the escape time i e 



l.hR/c. 
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Figure 1. Evolution of the particle distribution N(-y) correspond- 
ing to an injection of particles distributed in energy as a Gaussian, 
centered at 7 = 10 , for ti„j = O.lR/c. Labels indicate time after 
the beginning of the injection, in units of R/c (in the comoving 
frame). The escape time scale is t eS c = 1.5-R/c. The energy range 
has been divided according to the relevant time scales (top labels, 
see text). 



Figure 2. Evolution of the synchrotron self-Compton spectrum 
emitted by the particle distributions shown in Fig. 1, during a) the 
injection phase, and b) for t > ti„j. Labels indicate the time after 
the beginning of the injection, in units of R/c (in the comoving 
frame). These spectra do not take into account light crossing time 
effects. 



This is the value of t eac used in all simulations, unless oth- 
erwise noted. 

This case could represent a physical situation of a nar- 
row region of the jet invested by a perturbation which acti- 
vates the region. In this case R corresponds to the dimension 
perpendicular to the jet axis, while the width of the region 
could correspond to ~ t in jC <C R. Besides possibly repre- 
senting a real situation (even if highly idealized), this case 
is useful to clearly show the light crossing time effects. A 
more detailed discussion of the simulation of a region pow- 
ered by a shock will be presented in section 3.4. 

In order to better understand the different behaviours 
of the evolution of both particle and photons distributions, 
we can roughly divide the electron distribution and the emit- 
ted synchrotron spectrum into three main energy ranges 
depending on the characteristic involved time scales, with 
respect to R/c: 

(i) low energy (7 < 2 x 10 2 , v S yn ~ 10 11 Hz), in which the 
cooling time is at least one order of magnitude longer than 
R/c, so particle escape is the leading effect; 



(ii) medium energy 



lO^Hz, 



10 1J Hz) 



(2 x 10 2 < 7 < 2 x 10 4 , 
in which the electron cooling 



time is comparable to R/c, so that the cooling time is com- 
parable to the escape time; 

(iii) high energy (7 ^ 2 x 10 4 , v ayn > 10 15 Hz) in which 
tcooi <C R/c (by at least one order of magnitude), so that 
the particle distribution evolves more rapidly than R/c. 

In Fig. [l] we plot the electron distributions N("/) at 
different times, with high temporal resolution for the first 



5 - 



o 
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Figure 3. Evolution of the synchrotron self-Compton spectrum 
emitted by the particle distributions shown in Fig. 1, taking into 
account light crossing time effects. Note that the intensity is 
still rising at t = 0.5 R/c for frequencies under both the syn- 
chrotron and the Compton peaks, although the injection stops at 
t = 0.1 R/c. 
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instants, and lower temporal resolution for the rest of the 
evolution, when the cooling effects are much slower. 

The higher energy part of the distribution reaches the 
equilibrium state in a very short time, and does not vary 
during the injection time, except for a small decay due to the 
increasing radiaton field in the source, which increases the 
particle cooling. This is shown by the first three distributions 
in Fig. |l| labelled according to the time (in R/c units) after 
the beginning of the injection. 

Fig. ^| shows the synchrotron self-Compton spectra pro- 
duced by the distributions in Fig. [jj from t — 0.02 to 
t = WR/c without considering light travel time effects. The 
first instants (while the injection is active) are shown in Fig. 
|i]a, and the decay in Fig. gb. At frequencies above the syn- 
chrotron peak the cooling time is short, the equilibrium state 
is reached quickly and the flux is steady during the injection 
time,except for the very initial instants, not reproduced in 
Fig. Ha. Note that since electrons are injected only at very 
high energies their emission is first concentrated at high fre- 
quencies; only after some cooling time they can substantially 
emit at lower frequencies, simulating a newparticle injection 
at lower energies. This can be seen in Fig. Qb, where the flux 
at v < 10 13 Hz increases for t > tmj- 

Fig. ^ shows the spectra produced by considering the 
effect of t he l ight crossing time: as we already emphasized 
in section 2.3 this effect must be included when the N("y) 
distribution evolves on time scales shorter than R/c. In this 
case the observer will see different emissions coming from 
different regions of the source, as the source itself were not 
homogeneous. We stress that this is only due to the photon 
crossing time. In order to properly calculate the observed 
spectral evolution, we have to sum the correct contributions 
of the different source slices. The spectra shown in Fig. ^cor- 
respond to the first instants of the evolution. Note that the 
intensity of synchrotron radiation is now increasing during 
the injection time also at frequencies above the peak. This 
is mainly due to the photon crossing time: when the elec- 
tron cooling time is very short with respect to R/c (as it is 
for example for electrons emitting at the synchrotron peak 
at v ~ 10 17 Hz), the equilibrium state is reached quickly, 
and during the injection the observer receives photons com- 
ing from an increasing volume. After t = ti n j the flux re- 
mains steady until t = R/c, as can be easily seen also in the 
light curves (§ fc.l.l[ Fig. ^|). This effect is also shown by the 
Compton component. 

Due to the rapid cooling, electrons at the highest energy 
reach equilibrium in a short time, and during the injection 
time the distribution at these energies remains steady. After 
ti n j, it decays very rapidly. The corresponding synchrotron 
emission is then 'switched on' for a short time (during which 
it is constant) and then it is 'switched off'. At any given 
time (within R/c) the observer will then see a constant flux, 
produced by a single 'switched on' slice (which is 'running' 
across the source). This behaviour can be also seen in Fig. 
^ where the light curves at different frequencies are shown. 

3.1.1 Light curves 

We can compare the light curves obtained considering or 
neglecting light crossing time effects. This is shown in Fig. 

and Fig. ^p, respectively. We choose three frequencies on 
the synchrotron component as they are characteristic of the 
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Figure 4. Light curves of the specific intensity at different fre- 
quencies, corresponding to the spectral evolution of Fig. 2, to 
illustrate light crossing time effects, included in a) and ignored in 
b). Labels correspond to the logarithm of the frequency, as seen in 
the comoving frame. For clarity, the intensity of each light curve 
has been multiplied by different constants. 



three different energy ranges, and an additional frequency 
around the (initial) Compton peak. 

Fig. shows that at very low frequencies, (i.e. v ~ 10 10 
Hz), where t coo i 3> R/c, the rise of intensity is slow, because 
it is mainly controlled by the the light crossing time R/c: 
the observer will see the first slice, then the first and the 
second slices, etc.; when the last slice become visible, the 
first one is still emitting, because particle cooling is very 
slow. Since electrons are injected mainly at high energies, 
they can substantially emit at lower frequencies only after 
some cooling time scales: after t = R/c the intensity is still 
rising in some slices, causing the flux at these low frequencies 
to peak at t > R/c. 

The intensity at higher frequencies peaks instead at ear- 
lier times. At very high synchrotron frequencies [y ~ 10 16 Hz 
and above) we will have a plateau, because both t coo i and 
ti„j are much shorter than R/c: the electrons cool so fast 
that we will see a single slice 'running' across the source, as 
mentioned above. 

The shape of the decay phase is different for different 
frequencies: it is mainly controlled by t esc at very low fre- 
quencies (see the light curve at v — 10 10 Hz), while cooling 
dominates at higher frequencies. 

Note that the decay in the low frequency light curves in 
Fig. ^a are very similar to those in Fig. ^b, while the entire 
high frequency light curves and the initial instants of the 
low frequency ones cannot be reproduced without including 
the light crossing time effect. This happens because in these 
last two cases R/c is longer than t coo i or the injection time, 
respectively. 

In Fig. ^ we also show the u = 2 x 10 24 Hz light curve, 
which is close to the (initail) peak of the Compton emission 
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component. The shape of this curve is very similar to that 
of the v = 2x 10 16 Hz one, which is close to the synchrotron 
peak, because emission mainly comes from electrons of ap- 
proximately the same energy. 



R/c 



3.2 Power— law injection with ti„ 



R/c 



We now show the case of electrons distributed in energy as 
a power law Q(j) oc j~ p (p > 0), continuously injected for 
tinj ~ R/c. Input parameters for this case are: 7? = 10 16 
cm, Q(-f) oc 7" 1 ' 7 cm" 3 s" 1 with j min = 1, j max = 10 s , 



10" 



5 = 1 Gauss. 



We plot the time dependent spectra from t = R/c to 
t = 3R/c in Fig. Ha without considering the effect of the 
slices and in Fig. fib we summed the correct contribution of 
the different slices. 

Also in this case we roughly divide the synchrotron com- 
ponent into three energy ranges: 

(i) for v < 10 Hz the corresponding electron cooling 
time is much longer than R/c, and the decay of the flux (for 
t > ti„j) is mainly due to particle escape; 

(ii) for v > 10 Hz we have t coo i <C R/c: this is the part 
of the emitted spectrum in which we can see the greatest 
differences between the behaviours reported in Fig. [|a and 
|H|b; in other words, in this range of frequencies it is strictly 
necessary to take into account the different photon crossing 
times in order to correctly reproduce the evolution of the 
emitted spctrum; 

(iii) in the middle region we have t coo i ~ R/c. 

It is interesting to compare the flux behaviour, after 
the stop of the injection, in the spectral region in which 
tcooi -C R/c: in Fig. ^a (light crossing time effects ignored) 
we see that at the highest synchrotron frequencies the decay 
is very rapid, since the entire source closely follows the decay 
of the corresponding electron distribution. 

In the second case (Fig. |5|b, taking into account light 
crossing time effects) the behaviour is more complex, and 
after the end of the injection the emission falls showing three 
different phases: 

• t ^ R/c: we have a slow decay; we receive the emission 
from the distribution that starts to cool from the front slices, 
while the back slices are still 'switched on'. 

• t ~ 2R/c: when almost the entire source is (seen) 
cooled, the decay speeds up, because the only important 
contributions are given by the back slices, 

• t ~ 3R/c: the entire source will not contain highly ener- 
getic electrons any longer, and the emission at v > 10 15 Hz 
is dominated by the Compton component, produced by elec- 
trons of relatively small energy, whose cooling time is long. 
Correspondingly, at these frequencies the decay is again 
slow. 

This behaviour is illustrated by the light curves shown in 
Fig. | 

For the flux at the highest (Compton) frequencies, the 
behaviour is similar to the synchrotron flux, up to t ~ 2R/c, 
after which the flux disappears. 

It is important to note that just after the end of the 
injection, when high energy electrons are present and syn- 
chrotron photons are scattered at v > 10 24 Hz, we find sim- 
ilar variability factors for frequencies above the two peaks, 
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Figure 5. Evolution of the synchrotron sclf-Compton spectra 
calculated assuming to inject a power-law distribution, for ti n j = 
R/c. In a) light crossing time effects have been neglected, while 
in b) they have been taken into account. For clarity, only the 
decaying phase is shown. The shown spectra are separated in 
time by At = 0.25-R/c. Top labels divide the frequency range 
according to the relevant time scales (see text). 



contrary to the 'normal' behaviour which is characterized 
by a quadratic dependence of the Compton flux with re- 
spect to the synchrotron flux. This is due to the fact that 
the Compton scattering at these frequencies is in the Klein 
Nishina regime: high energy electrons efficiently scatter seed 
photons of frequencies hv/m e c 2 < 1/7 (see equation 16). 
For the highest energy electrons the available synchrotron 
photons for scattering are less, reducing the variability am- 
plitude (see e.g. Ghisellini, Maraschi & Dondi 1996). 



3.2.1 Light curves 

We report in Fig. ^a and |^b the light curves at four differ- 
ent frequencies considering or not light crossing time effects. 
We can see that taking them into account we can have sym- 
metric light curves, with equal rising and decaying phases, 
at frequencies for which t coo i <C R/c. The cooling time of 
the electrons emitting at high frequencies (10 15 — 10 17 Hz) 
is much shorter than R/c, and the electron distribution re- 
mains steady for the entire injection time. This implies that 
for t < R/c the observer sees an increase of the flux due 
to the increasing 'switched on' volume of the source. Af- 
ter t = R/c, the front slices 'switch off', and the total flux 
decreases. This corresponds to symmetric light curves with- 
out plateaux. At lower frequencies, at which t coo i > R/c, 
the decay is slower than the rise, and the light curves are 
asymmetric. 

Another remarkable effect is the existence of time-lags 
among the light curves at different frequencies: it is easy to 
note that the emission peaks are not reached at the same 
time and the flux at the highest frequencies appear to lead. 
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Figure 6. Light curves of the specific intensity at different fre- 
quencies, corresponding to the spectral evolution of Fig. 5, to 
illustrate light crossing time effects, included in a) and ignored in 
b). Labels correspond to the logarithm of the frequency, as seen in 
the comoving frame. For clarity, the intensity of each light curve 
has been multiplied by different constants. 



All these behaviours are not observable if we do not 
take into account the light crossing time: in this case we 
see the formation of a plateau at high frequencies, where 
the equilibrium state is reached first. Furthermore we can- 
not observe time lags between different frequencies, because 
when the injection stops, the intensity begins to fall at all 
frequencies, although with different decay shapes (steeper at 
higher frequencies). 

In Fig. we can observe the changes in the slope of the 
decay of the v = 10 16 3 Hz light curve, according to what 
discussed in Sect. 3.2. 



3.3 Time delays 

We have already noticed the existence of time-lags among 
light curves at different frequencies, that are clearly visible 
in Fig. ^|a: this is a consequence of taking into account the 
light crossing time effects, although time delays are mainly 
caused by the different cooling times of electrons emitting 
at different frequencies. 

Assume in fact that tinj = R/c: at t = tinj the observer 
sees the stop of the particle injection in the 'front' of the 
source, while the back of the source is still 'switched off'; in 
other words, at this time the observer is receiving the spectra 
emitted by a population of electrons at t evo i = R/c (from 
the front of the source) and at t evo i = (from the back of the 
source) . At later times in the front of the source the electron 
distribution starts to decay (with a corresponding decay of 
the emitted flux), while in the back the flux is still rising 
until the particle distribution reaches, after some t coo i, the 
equilibrium state. This combination of slices in which the 
emitted flux is increasing and slices in which the decaying 



phase is already started, determines the position of the peak 
on each light curve. Notice that if the electron cooling time 
is very short (for example in the case of the v — 10 16 Hz in 
Fig. H) , the equilibrium state is reached in a very short time 
with respect to R/c: after the injection stops, each slice is 
quickly turned off, and the peak of the light curve is reached 
soon after t = ti n j = R/c. In light curves corresponding to 
lower energy electrons the peak will occour later, because 
the equilibrium state is reached later: after the stop of the 
injection the flux is slowly decreasing in the front of the 
source, while it is still increasing in the back. The longer 
the particle cooling time, the longer the time that the light 
curve needs to reach the maximum. 

We emphasize that the different cooling times of the 
electrons responsible for the emission at different frequen- 
cies introduce more visible time-lags if the light crossing 
time effects are included, as can be seen comparing the light 
curves for the high synchrotron frequencies shown in Fig. 6a 
and Fig. 6b. 

3.4 Simulation of a shock 

We assume that a shock of longitudinal dimension R and 
width r <C R runs along a region of the jet of same dimension 
R (perpendicular to the jet axis). Particle flows along the 
jet with a velocity (3c. In the frame comoving with this flow 
the shock moves with a velocity j3' s c. As in the previous 
sections, we assu me that the observer is located at an angle 
1/r = yl fl 2 from the jet axis, such that the viewing 
angle in the comoving frame is 90°. We assume that the 
shock is active for a time t s , as measured in the comoving 
frame. At a given location, the shock accelerates particles 
for a time t in j = r/(f3' g c). 

To show how to correctly sum the different contribu- 
tions of different emitting regions, we sketch in Fig. (?] an il- 
lustrative example with a small number of different regions. 
In this case, each slice parallel to the jet axis is characterized 
by different electron distributions: 'young' electrons close to 
the shock front, and older ones far from it. It is therefore 
necessary to subdivide each slice into an appropriate num- 
ber of 'cells' (see Fig. Q). The velocity of the shock in the 
frame comoving with the particle flow determines the size 
of the cells in the direction parallel to the jet axis, and then 
the total emitting volume. 

Fig. ^a shows the light curves obtained in the case of a 
shock of width r = tinjP' s c, with tinj = O.lR/c active for a 
time ts = R/c, with all other parameters equal to the case 
of §EO. 



Comparing Fig. 8a with Fig. 8b one can notice that 
in the shock case we obtain further delays between different 
frequencies, due to the different way of summing the spectra 
produced at different timesteps. 

From these simulations we can conclude that there are 
no qualitative differencies between the light curves in the 
'homogeneous' case with ti n j = R/c, except for the occur- 
rence of further time delays between light curves at different 
frequencies. 

3.5 Multiple rapid injections 

In Fig. ^a we show the light curves resulting from discontin- 
uously injecting a gaussian distribution of particles, centered 
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Figure 7. Schematic illustration of the composition of spectra 
arriving at the observer at the same time, but leaving each shell 
at a different time, hence produced by electrons of a different 
age. The observer is assumed to be located at 90° to the shock 
velocity (bottom of the figure) in the frame comoving with the 
particle flow. In the left panel the shock has been active for a 
time R/f3' s c since the beginning of the injection. The number in 
each cell indicates the "age" of the electrons. Note that above the 
"diagonal" line (cells without number), photons had not enough 
time to reach the observer. The right panel shows how to sum 
the contributions of the cells at the time equal to 2R/(3' s c, in the 
case that the shocks stops to inject particles after t s = R/c. At 
later times each cell will contribute to the emission with an older 
distribution, and correspondingly the number characterizing each 
shell will increase by one unit at each timestep At. 



at 7 = 10 s . In the shown case, we have assumed that there 
are 17 injection phases, lasting from tinj = O.lR/c to tmj = 
R/c, separated by different times, and that t esc — 8R/c. The 
electron injection rates are equal for all cases. 

It can be seen that these injection phases can be re- 
solved at the highest synchrotron frequencies, while they 
merge at lower frequencies, where the memory of the differ- 
ent injection phases is lost. At energies at which t coo i > R/c, 
the electron distribution builds up, since electrons cool in a 
time scale longer than the time separation between different 
injection phases. This introduces both a smoothing of the 
corresponding light curve and a time delay with respect to 
higher frequency light curves. 

The series of short injection phases separated by short 
times makes the synchrotron flux at the highest frequencies 
to flicker. The superposition of the different distributions 
corresponding to each injection phase prevents the forma- 
tion of the 'plateau', characteristic of a single short time 
injection. This also causes the variability time scales to be 
shorter than R/c. 

In Fig. ^|b we show the time evolution of three spectral 
indices (an. 5-12,6, a12.s-14.fi and a14.5-i6.fi, connecting the 
flux at log v — 11.5 and log v — 12.5, the flux at log v = 12.5 
and log v = 14.5 and the flux at log v — 14.5 and log v — 16, 
respectively, where frequencies are measured in Hz). The 
spectral indices can vary (even by a large amount) during the 
major injection phases, while they remain within a narrow 
range during the flickering phases. 
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Figure 8. a) Light curves at different frequencies (as seen in 
the comoving frame) corresponding to a simulation of a shock of 
width r = t{ n j(3' s c, t{ n j = 0.1-R/c, injecting energetic particles 
for a time t s = R/c. All other parameters are the same as in 
the simulations shown in Fig. 6. In panel b) we report the the 
light curves already shown in Fig. 6a in the case ti n j = R/c. For 
clarity, the intensity of each light curve has been multiplied by 
different constants. 



4 DISCUSSION 

We summarize here the main characteristic of the light 
curves for different values of tinj, t CO oi and t esc : 

• t in j < R/c 

- when tc 00 i S> R/c we see a slow rise of intensity, 
mainly controlled by the light crossing time R/c: the 
maximum of intensity is reached for t > R/c, and de- 
pends on t coo i and t esc . At small frequencies, where 
tcooi S> tesc, the decay of the emission is controlled by 
t eac , while at frequencies where t coo i < t esc the leading 
effect is cooling; 

- when t coo i — R/c the rising shape of the light curve 
is similar to the above case. The maximum of intensity 
is now reached at t ~ R/c, and the shape of the decay 
is controlled by t coo i', in the special case in which R/c ~ 
some tcooi we have a symmetric light curve; 

if tcooi ti n j we have a fast rise of intensity (t r ise — 
ti„j) and then a plateau (or a very slow rise in the shock 
case) until t ~ R/c. For t > R/c the decay is as fast as 
the initial rise. 



* tin 



R/c 



- if tcooi ^ R/c we see a slow rise, the maximum is at 
t w 1.5-R/c, and the decay is due to escape if t coo i 3> t esc ', 

- if tcooi ~ tinj ~ R/c the shapes of the rising and de- 
caying phases are very similar, resulting in a symmetric 
light curve; 

- where tcooi R/c both the rise and the decay are 
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Figure 9. Light curves of the specific intensity at different fre- 
quencies (as seen in the comoving frame) and spectral indices 
resulting from multiple injection phases. For clarity, the intensity 
of each light curve has been multiplied by different constants. 



controlled by the light crossing time, and we have sym- 
metric light curves with time lags depending on the dif- 
ferent cooling times; 

• if ti n j is much longer than R/c and t coo i the light curve 
will have a plateau, because the long injection time allows 
the entire source to reach equilibrium. 

The main difference between injecting electrons at high 
energies ('monoenergetic' or gaussian injection) and inject- 
ing a power-law (with p > 0) distribution is that in the first 
case the emission will be concentrated first at high frequen- 
cies, and only after some t coo i electrons can substantially 
emit at lower frequencies. This simulates a sort of a 'new' 
injection at low energies, that continues for t > tinj, and 
produces a time delay between the peak of the emission at 
different frequencies. 

If the emitting plasma has bulk motion, we must take 
into account the effects of beaming. If T is the bulk Lorentz 
factor, 6 the viewing angle and S — [F(l — /3cos6)] _1 the 
beaming factor, the observed intensity is I(v) = S 3 I'(i/) 
and t = t' /S, where I'{v') and t! are the comoving intensity 
and comoving time scales, respectively. 

Observations of variability in BL Lac objects suggest 
the simultaneous presence of different variability time scales 
on each source (e.g. Ghisellini et al. 1997, Massaro et al. 
1996), suggesting that the variations can be originated by 
more than one component: a larger one can take into account 
of the longer time scale variations, and smaller ones can 
originate the rapid flares. Flares are observed at almost all 
wavelengths (depending on each source characteristics), and 



can have different durations (from hours to days). The oc- 
currence of outbursts with symmetric shape (similar rise and 
fall) have been recently reported for well monitored sources, 
especially in the optical band (Massaro et al. 1996, Ghis- 
ellini et al. 1997) and in the X-ray band (Urry et al. 1997). 
According to our results, this behaviour can be originated 
only in two cases: 

(i) t in j "C R/C ~ tcooi 

(ii) ti n j ~ R/c and t coo i < R/c. 

In case (i) symmetric light curves are present only within a 
very small range of frequencies (t coo i ~ R/c), while in case 
(ii) quasi-symmetric light curves can occur at all frequencies 
corresponding to particle cooling time scales shorter than 
R/c. Furthermore the second case can be interpreted as a 
result of a shock lasting for a time tinj ~ R/c. 

Note that the underlying component has not been in- 
cluded in the previous simulations (§§3.1 and 3.2), which 
only describe the evolution of the 'flaring' component. To 
correctly reproduce the blazars variability behaviour, it is 
therefore necessary to take into account the presence of a 
quasi-stationary component, diluting the variability of the 
flaring component. 

Finally, note that the simulations in §3.4 of a shock 
active for a time R/c are similar to the corresponding 'ho- 
mogenous' case, with slightly longer time delays between 
light curves at different frequencies. 



5 APPLICATION TO MKN 421 

Mkn 421 is one of the nearest BL Lac Lac objects (z = 0.03) 
and it is classified as XBL (X-ray selected BL Lac) (Giommi 
et al. 1990, Hewitt & Burbidge 1993). In the overall spectral 
energy distribution [v — vF(y)\ the X-ray emission smoothly 
connects to the optical and UV. It is generally believed that 
this component is due to synchrotron emission. Mkn 421 
is a faint EGRET source with a flat GeV spectral index 
(Sreekumar et al. 1996) and it was the first blazar detected 
at TeV energies (Punch et al. 1992). In the optical, UV, X- 
ray and TeV bands it shows strong and rapid flux variability 
(Buckley et al. 1996). 

In May 1994 the ASCA satellite revealed an X-ray flare 
(Takahashi et al. 1996) during an high state of TeV emission 
(Macomb et al. 1995). Observations report an increase of a 
factor ~ 2 of the 2-10 keV flux, with a doubling time scale of 
~ 12 hours, while much less amplitude variability is present 
in the IR, optical, UV and GeV (EGRET) bands. Takahashi 
et al. (1996) found a time-lag between hard X-rays (2 — 7.5 
keV) and soft X-rays (0.5 — 1 keV) of ~ 1 hour: the hard 
X-rays lead the soft X-rays. They interpret this as due to 
synchrotron cooling. With the numerical code described in 
the previous sections we can qualitatively reproduce this be- 
haviour, first by fitting the high and the low states (Macomb 
et al. 1995) as equilibrium SSC spectra, and then by gen- 
erating simulated light curves at all frequencies. In Fig. |l(j 
we report the simultaneous spectra taken during the high 
and low states (data from Macomb et al. 1995, 1996, using 
Ho — 50 Km s _1 Mpc -1 , go = 0.5). The parameters used 
for the fits are reported in the figure caption. 

We assumed that the rapid variability revealed by 
ASCA during the highly active X-ray/TeV state is due 
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Figure 10. High and low states of Mkn 421 during the 1994 mul- 
tiwavelength campaign (Macomb et al. 1995); continuous lines 
are equilibrium SSC models. We found satisfactory fits with the 
following parameters: R, = 4.8 X 10 16 cm, B = 0.08 Gauss, in- 
jection law Q(7) OC 7 -1 ' 5 exp( — ~t/~jmax) with "/ min = 1, 7mm = 
1.3 X 10 5 , t e3c = 4R/c, injected compactness £{ n j = 9 X 10 -5 , 
beaming factor 5 = 15.5 (low state); R = 4.8 X 10 16 cm, 
B = 0.04 Gauss, (3(7) oc 7 -1 ' 7 exp(-7/7 maa; ) with 7 min = 1, 
■y max = 8.5 x 10 5 , £ inj = 2.5 x 10" 4 , t esc = 3R/c, 5 = 15.5 (high 
state). 



to the sum of a rapidly evolving component and a quasi- 
constant one, corresponding to the high state fit. We add to 
this constant emission a variable component, for which we 
found the following parameters: 

R = 1.5 X 10 16 cm, B = 0.13 Gauss, S = 15.5, £ inj = 
1.5 x 1CT 3 , Q(7) oc 7 1 ' 4 exp(~-y/-y max ) between -y min = 10 3 
and 'jmax = 8.5 x 10 5 . We perform the simulation in the 
shock case with the following parameters: r s = O.lR/c, t s = 
R/c and f3' s ~ 1. 

We plot the obtained light curves at four different fre- 
quencies in Fig. jnj. Note that the flux is still increasing at all 
reported frequencies for t > t s (l + z)/S, which corresponds 
to t ~ 9 hours: this is due to the different evolving state of 
the electrons in the different cells. 

Because of the photons travel time, when the injection 
stops (at t s = R/c) the observer still sees the emission from 
the back of the source produced by the increasing parti- 
cle distribution. Therefore the maximum flux is observed at 
a time greater than R/c, by an amount depending on the 
cooling time: the flux at the highest frequencies peaks first, 
followed by the flux at lower frequencies. As a consequence, 
the flux at the synchrotron peak (X-ray band) leads the flux 
at lower frequencies, and time-lags can be observed even be- 
tween close frequencies within the same band. 

Another consequence of our modelling is the prediction 
of time-lags between the synchrotron and the self-Compton 
fluxes. While the X-ray and the TeV emission are proba- 
bly produced by the same electrons, the most effective seed 



Figure 11. Ligh curves at different energies, reproducing the 
rapid flare shown by MKN 421 on May 16, 1994. In a) we plotted 
the flux at different frequencies normalized to the steady compo- 
nent; in b) the spectral index connecting the 4.3 keV and the 0.9 
keV fluxes is shown. Times and frequencies take into account of 
the beaming factor <5. 

photons for the TeV emission are in the (observed) optical- 
UV band, due to the Klein Nishina decline of the scattering 
cross section. In the case discussed here, the electron pop- 
ulation producing the optical-UV emission builds up, as a 
consequence of particle cooling, some time after the begin- 
ning of the injection. Optical-UV and TeV emissions should 
therefore be nearly simultaneous, lagging the X-ray emis- 
sion. A delay between X-ray and TeV emission can also be 
introduced by taking into account the different light travel 
times of the seed photons to be scattered at high energies, 
as mentioned in §2.3, but to be more quantitative we must 
await a more detailed numerical treatment. 

A time lag between X-ray and TeV fluxes should exist 
only if the seed photons for the scattering are produced lo- 
cally. In the alternative case of seed photons coming from 
external regions (by, e.g. a dusty IR torus), we expect to 
observe a simultaneous rise in the X-ray and TeV fluxes of 
the same factor. Note that the 'external' radiation could also 
come from adiacent regions of the flaring component, i.e. the 
ones responsible for the 'diluting' component. Also in this 
case the X-ray and the TeV fluxes, being produced by the 
same electrons, should vary linearly and simultaneously. 



6 CONCLUSIONS 

By summing the contribution of different slices of the source, 
the observer receives photons produced by particle distri- 
butions of different ages, a situation which resembles the 
one occuring in an inhomogeneous source. This effect is im- 
portant whenever the particle distribution evolves on time 
scales shorter than R/c. Despite the fact that this effect 
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tends to smooth out fast variations, time delays between 
light curves at different frequencies are observable, as illus- 
trated for the hard/soft X-ray emission of Mkn 421. Our 
results are particularly important for the observed fast vari- 
ability of blazars, where the variability time scales indicate 
extremely compact emission regions, thus large magnetic 
and radiation energy densities, and consequently cooling 
time scales shorter than the light crossing time R/c. This 
is also a necessary condition for having symmetrical flares 
(equal rising and decay time scales), often observed in the 
optical and in the X-ray band. Note that the case of 3C 
279, one of the best studied blazars which showed asym- 
metric light curves at high energies, is probably different, 
because the high energy emission can be inverse Compton 
scattering off seed photons produced externally to the jet, 
and the resulting light curve can be connected to the geo- 
metrical distribution of these seed photons: for example, if 
an active blob passes through the broad line region (BLR), 
there will be enhanced inverse Compton emission as long as 
the blob is inside the BLR, and a sudden drop when the 
blob moves outside the BLR. 

Within the synchrotron self-Compton scenario, the 
symmetrical behaviour of the light curves near or above the 
peak of the synchrotron emission can be contrasted with 
the light curves at lower frequencies, characterized by cool- 
ing (and escape) time scales longer than R/c. In this case 
the light curves should be asymmetric, with a rising phase 
lasting £ R/c (if U n j ~ R/c), and a decaying phase lasting 
for t coo i (or t esc , if shorter than t coo l)- Sources with different 
magnetic field and/or different compactnesses should behave 
differently, since for large compactness (or larger magnetic 
field) sources, the cooling time can be shorter than R/c for 
a wider range of electron energies, resulting in symmetric 
light curves for a wider range of frequencies. 

We stress that symmetric light curves without plateaux 
strongly constrain the injection and the cooling time scales: 
if both these times are much longer than the the light cross- 
ing time R/c the light curve is not symmetric; if both time 
scales are much shorter than R/c the curve is symmetric, but 
a plateau forms, while if the injection time is shorter than 
R/c and the cooling time is of the same order, a symmet- 
ric light curve is possible (without plateau), but only at one 
specific frequency. Therefore symmetric light curves without 
plateau at more than one frequency are possible only if the 
the injection last for R/c and the cooling time is shorter 
than R/c. 

A more complex behaviour is possible if the injection 
of relativistic electrons is impulsive and repeated several 
times within one light crossing time R/c. In this case the 
synchrotron flux at the largest frequencies responds to the 
different injection phases, and the repeated injections can 
make the plateaux in their light curve to disappear. Varia- 
tions of relatively small amplitude are possible in very short 
times, which would lead to calculate variability time scales 
shorter than R/c. At lower frequencies, where cooling times 
are longer, the electron distributions corresponding to dif- 
ferent injections can build up and the memory of the in- 
dividual injection phases is lost. At these frequencies the 
light curves are smoother, and they can have apparent de- 
lays (with respect to higher frequency light curves) produced 
by this 'build up' effect. 
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